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Abstract. We use nonequilibrium renormalization group (RG) techniques to analyze 
the thermalization process in quantum field theory, and by extension reheating after 
inflation. Even if at a high scale A the theory is described by a non-dissipative Aip 4 
theory, the RG running induces nontrivial noise and dissipation. For long wavelength, 
slowly varying field configurations, the noise and dissipation are white and ohmic, 
■ respectively. The theory will then tend to thermalize to an effective temperature given 

by the fluctuation-dissipation theorem. 
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^ . The goal of this paper is to show how nonequilibrium renormalization group (RG) 

techniques may be applied to study the thermalization process in quantum field theory 
and, by extension, the turbulent reheating period in inflation. 

The issue of thermalization in quantum field theory has received renewed attention 
in later years, motivated by applications to cosmology and to relativistic heavy ion 
collisions, as well as by theoretical progress through large scale numerical experiments 
PQ. However, there is a lack for analytical, model robust methods capable of yielding 
predictions for such observables as the final temperature and thermalization time scales. 
This need is particularly felt in applications to cosmology, since reheating is more 
likely than not a complex phenomena involving several nonlinear fields and the evolving 
background geometry [2]. 

During inflation, the dominant form of matter in the Universe is a condensate 
(the inflaton) which evolves rolling down the slope of its effective potential. When the 
inflaton nears the bottom of the potential well, it begins to oscillate and transfers its 
energy to ordinary matter (then in its vacuum state). We call this process reheating. 
Reheating proceeds through several stages [3] correlated with the different stages of 
the thermalization process, namely preheating, inflaton fragmentation and turbulent 
thermalization. Generally speaking, the early phases produce an spectrum with high 
occupation numbers in a narrow set of modes. Turbulent thermalization concerns the 
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spread of the spectrum over the full momentum space and the final achievement of a 
thermal shape. 

At early times occupation numbers are high and the process may be described in 
terms of classical wave turbulence [3j. As the spectrum spreads occupation numbers 
fall and the classical approximation breaks down. The challenge for us is to provide a 
quantum description of turbulent reheating. 

For demonstration processes, we shall only discuss quantum turbulent thermaliza- 
tion in a nonlinear scalar field theory in 3 + 1 flat space-time. 

The basic idea is the same as in Kolmogorov - Heisenberg turbulence theory: a mode 
of the field with wave number k lives in the environment provided by all modes with wave 
number k! > k. Since the physical mechanism for damping in the long wavelength sector 
is the interaction with shorter wavelength modes, it is natural to understand damping 
as a feature of the effective theory where the shorter modes have been coarse-grained 
away [5j[6]. Since this operation will leave the long wavelength modes in a mixed state, 
the natural description of the relevant sector is in terms of a density matrix [7], and 
the natural action functional encoding the effective dynamics is the Feynman- Vernon 
Influence Functional (IF) [HI [9] . Now suppose we are given the IF when all modes k > A 
have been coarse-grained away, and we wish to further coarse grain the modes in the 
range A > k > k . We split the desired range into shells of infinitesimal thickness 5s, 
and integrate out each shell retaining only terms of order 5s. Adding a change of units 
and a rescaling of the fields after each integration, we transform the shell coarse-graining 
into a RG flow in the space of influence functionals [1QJ. Because we shall not assume 
equilibrium conditions, this may be called the nonequilibrium RG. 

Here A is not meant as an ultraviolet cutoff to be removed eventually, but rather 
as the "hard" scale at which the microscopic theory is well understood and radiative 
corrections are perturbative. Our goal is to investigate physics at "soft" scales k <C A. 
This issue has been studied in the context of hot nonabelian plasmas, where the 
emergence of dissipation and noise has been demonstrated in different ways [TTl [T2l [T3] . 
We wish to emphasize the nonequilibrium aspects of the problem, as well as to put those 
findings on a more general base by adopting the RG approach. 

It is important to stress two basic differences between the nonequilibrium and 
equilibrium RG [Hj. The IF may be regarded as an action for a theory defined on 
a "closed time path" (CTP) composed of a first branch (going from the initial time 
t = to a later time t = T when the relevant observations will be performed -that 
is why we need the density matrix at T) and a second branch returning from T to 
[T5| [TBI [TTl [TH] . Thus each physical degree of freedom on the first branch acquires a twin 
on the second branch -we say the number of degrees of freedom is doubled. The IF is not 
just a combination of the usual actions for each branch, but also admits direct couplings 
across the branches. The damping constant k and the noise constant v are associated 
to two of those "mixed" terms. Therefore, the structure of the IF (from now on, CTP 
action, to emphasize this feature) is much more complex than the usual Euclidean or 
"IN-OUT" action. 
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The second fundamental difference is the presence of the parameter T itself. In 
nonequilibrium evolution, it is important to specify the time scale over which we shall 
observe the system. The CTP action contains this physical time scale T. From the 
point of view of the RG, this adds one more dimensional parameter to the theory, 
much as an external field in the Ising model. Physically, because time integrations are 
restricted to the interval [0, T], energy conservation does not hold at each vertex. This 
is of paramount importance regarding damping. 

The RG for the CTP effective action (obtained by taking the limit T — > oo) was 
studied by Dalvit and Mazzitelli [T9J [20] ; see also refs. [6] and [7] . Unlike those works, 
we focus on the dissipation and noise features of the effective dynamics, rather than in 
the running of the effective potential. 

In formulating a nonequilibrium RG, we must deal with the fact that the CTP 
action may have an arbitrary functional dependence on the fields and be nonlocal both 
in time and space. In principle, one can define an exact RG transformation [T9] , where 
all three functional dependencies are left open. However, the resulting formalism is too 
complex to be of practical use. Fortunately, the special properties of the application to 
thermalization allow for a substantial simplifications, such as working in three spatial 
dimensions. 

The full RG equations for this theory are given in [21]. Here we shall only 
highlight those aspects of the calculation most relevant to the application to turbulent 
thermalization. 

Let us call ^ the field variable in the first (resp. second) branch of the CTP. To 
write down the CTP action, it is best to introduce average and difference variables 

(f) = <p+-<p-, (1) 

<p = <p + + <p~. (2) 

In terms of these variables, a generic CTP action may be written as 

'S'CTP = Sq + S\ + ^her, (3) 

where Sq is the CTP action functional for a free massless field theory 



1 rT 



S [<p,<p] = -J Q dt Jd a k [0(k,t) <p(-k,t)-k 2 0(k,%(-k,t) 
S\ accounts for a A<^ 4 -type self interaction 
c u i A f T m [d d k 1 ...d% d ( 4 \ 



(4) 



48 Jo J (2n) d V/ .., ; 

x 0(k 1; *)p(k 2 , £Mk 3 , %(k 4 , t) + 0(k 1; t)0(k 2 , t)0(k 3 , %(k 4 , t)] , (5) 

and S'other includes all other possible terms. Momentum integrals are bounded by k = A, 
and d = 3. We shall assume that the initial condition for the RG flow is S'other = at 
the hard scale A, so that if it appears at soft scales, it is as a consequence of the RG 
running itself. Note that this is true, in particular, for the noise and dissipation terms. 
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To define the nonequilibrium RG we also need to specify the state of the field at 
the initial time t = 0. For simplicity, we shall assume this is the vacuum state for the 
free action Sq. Observe that this is a nonequilibrium state for the interacting theory. 

The value A of the coupling constant A at the hard scale A may be used as the 
small parameter in a perturbative expansion of the RG equation. To order Aq, the RG 
equation for the quartic coupling decouples, and can be solved by itself. The result is 
that at soft scales k, A is both scale and T dependent. There is no RG running if T = 0, 
while the usual textbook result is obtained as T — ► oo [22J. For all values of T, A is 
driven to zero as k — > [21]. Thus it is consistent to assume that A is uniformly small 
in the relevant scale range. 

In particular, in order to compute the RG equations to order Aq, it will be enough 
to use in the Feynman graphs the zeroth order propagators, which are those of the 
massless free theory. The only exception is in computing the effective mass, but this 
calculation is decoupled from the noise and dissipation terms to order Aq. Observe that 
it is at the same time a huge simplification and a strong limitation concerning the range 
of application of our results, as we expect substantial shifts in the propagators when T 
approaches the relaxation time of the theory. 

Because of the nonzero initial value of A, other couplings will appear as a result of 
the RG running. To order Aq, it is enough to consider quadratic, quartic and six-point 
terms in the action. All these terms feed back into each other, so they must be taken 
self-consistently. If we understand thermalization in the usual sense that propagators 
become approximately thermal [23], however, it is enough to focus on the quadratic 
terms, 

Mother -> S 2 [(f),(p) = f dti f dt 2 f d d k \v 21 (k;t x ,t 2 ) (j)(k,t 1 )<p(-k,t 2 ) 
Jo Jo J 1 

+ ii>22(fc;ti,*2)0(Mi)0(-k,*2)]. (6) 

In principle, the induced quadratic terms will be oscillatory functions of At 12 . 
However, we are interested mostly in the dynamics of slowly varying field configurations 
which are insensitive to high frequencies. To focus on the slow dynamics, we may project 
out the mass, dissipation and noise terms on which the oscillations are mounted. 

To this end, we introduce two projectors. Given a function of two times v(k; t\, t 2 ), 
we define 

Pv(k;t 1 ,t 2 )=Vv(k)5(t 1 -t 2 ), (7) 
and, if v(k; ti, t 2 ) = for t 2 > t\, 

Qv(k;t u t 2 ) = Qv(k) 



2^ + 5(t 2 )-5(0)j 



(8) 



where 

T Jo Jo 



Vv{k) = 1 f dt x ( T dt 2 v(k; h, t 2 ), (9) 
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and 

Qv(k) = )- ( T dt, ( T dt 2 v(k; h, t 2 ) {t 2 - h). (10) 
1 Jo Jo 

It is easy to verify that P 2 = P, Q 2 = Q, and that QP = PQ = 0. This proves that 
the decomposition 

v(k; h, t 2 ) = Pv(k; t h t 2 ) + Qv(k; t u t 2 ) + Av(k; t h t 2 ) (11) 

is unique. When this decomposition is applied to v 2 \ in equation ([6]), we get 

F dt x [ T dt 2 [ d d kv 21 (k;t u t 2 ) 0(k,*i)p(-k,t 2 ) = 
Jo Jo J 

f dtx J d d k [v (k) 0(k, hM-k, ti) + v x (k) 0(k, tjtpi-k, h)] 

f T dt x f T dt 2 f d d k Av 21 (k;t 1 ,t 2 ) 0(MiM-M 2 ), (12) 
Jo Jo J 
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where 



and 



^ „x / x , ^o(O) k 2 d 2 v (0) , , 

«b(*) = = «b(0) + + y.^W 1 + ■ ■ ■ . ( 13 ) 



^(A) = Qu 21 (A0. (14) 

The linear term in equation ( 1T31) vanishes from symmetry, and the appearance of 
the quadratic term is prevented by performing a field rescaling as part of the RG 
transformation (thus the field acquires an anomalous dimension). Neglecting the last 
term in equation f|T2l . the net effect for the long wavelength modes is to induce a mass 
term m 2 = — 2t>o(0), and a damping constant n(k) = —v\(k). The noise kernel is 
obtained in a similar way from the imaginary part of the CTP action, v 22 in equation 
©. 

After these considerations, the relevant CTP action for long wavelength, slowly 
varying configurations reduces to 

S C tp[0, if] = [dt J d d k [±0(k, t) (p(-k, t) - i<Kk, t) (k 2 + m 2 ) (fi-k, t) 

- n(k) 0(k, tM-k, t) + l -v{k) 0(k, t)0(-k, *)] . (15) 

It is shown in Feynman and Hibbs [9] that this CTP action describes a field subject 
to ohmic dissipation with damping constant k and a stochastic source j (t) with white 
noise self-correlation (j (t) j (t')) = v 5 (t — t') . The relationship of the propagators of 
the original theory to those obtained from this CTP action is discussed in [21]. For 
present purposes, it is enough to observe that this system thermalizes to an effective 
temperature given by the fluctuation-dissipation theorem [8j [9j [23, [26] 

T« = £, (16) 
with a thermalization time 

r = -. (17) 

K 



Nonequilibrium renormalization group 



6 




Figure 1. The dissipation constant k (clashed), the noise v (dotted), both measured 
in units of (A/967T 2 ) 2 , and the effective temperature T e s as functions of the scale k for 
a fixed value of T = 0.5. 



For fc<A we obtain the approximate expressions [21] 



and 

v(k,T), 
+ 4 InT 




T) ~ 4 (9^) [? - 2T ' " 8 cos T + cos ( 2T )] ln ( A /*) 



135 + 4 [ 7b - 34 In 2 - 7 In 3 + 7 Ci(3T)l - 3 cos(4T) 



- 12 



I 6 8tt + 3T - 16 Si(2T) ] sin T + 4 [7tt + 6T - 14 Si(4T)j sin(2T) 
19 - 4T 21 



Ci(T) + 8 7 - 6T 2 Ci(2T) + 4 



1 -4T' 



Ci(4T) - 8 cos(3T) 



- 24[7-8Ci(2T)]cosT + 4 [ll + 2Ci(T) -2Ci(2T) - 14Ci(4T)] cos(2T) 

+ 4T[18tt - 54 Si(T) + 28 Si(2T) - 6 Si(3T) - 2 Si(4T) + sin(4T)]|, (19) 

where Si and Ci are the sin and the cos integral functions. The thermalization time and 
the effective temperature go as ln(A/A;) _1 when k — > 0. Observe that the asymptotic 
formula for k is not positive definite. This suggests that the fundamental damping 
mechanism is Landau damping of long wavelength waves through interaction with hard 
quanta [27J. In any case, we expect our approximations to break down before we reach 
the point k — 0. 

In figure 1 we show k, u, and T e ff as functions of the scale k for fixed T . In figure 
2 we show k, u, and T e g as functions of the observation time T for a fixed k. We have 
chosen units in such a way that A = 1. The expressions for k and v are given in |21j . 

The most important result of this paper is that RG running alone describes the onset 
of noise and dissipation in the long wavelength modes, even if these are not assumed to be 
present in the underlying microscopic theory. These two elements provide the sufficient 
conditions for thermalization as described by Schwinger [TJ)J [TBI El EH] • Therefore the 
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Figure 2. The dissipation constant k (dashed), the noise v (dotted), both measured 
in units of (A/967T 2 ) 2 , and the effective temperature T cff (solid), as functions of T for 
a fixed value of k = 1/60. Each quantity is normalized with respect to its maximum 
value in the displayed interval(w 215 for k, 26 for v and 10 for T e g). 

model succeeds in describing the onset of the thermalization process, and yields simple 
estimates of both the final temperature and the thermalization times at different scales. 
We must caution however that in this form these estimates are not fully reliable, as they 
involve pushing the theory to the range AT > 1. This lies beyond the range of validity 
of our approximations concerning T. To extend the nonequilibrium RG to a larger T 
range a fully self-consistent approach is necessary [28j [23 EQ] . 

Because of this limitation, we do not claim to have solved the problem, but to 
have shown a framework for a solution. We need a self-consistent approach to the hard 
loops to be able to extend further the T range. Still, the power of the nonequilibrium 
RG allows us to extend to the fully quantum regime the insights gained from wave 
turbulence in the classical stage of evolution. 
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